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Abstract 

Prompted by recent experimental developments, a theory of surface scattering of fast atoms at 
grazing incidence is developed. The theory gives rise to a quantum mechanical limit for ordered 
surfaces that describes coherent diffraction peaks whose thermal attenuation is governed by a 
Debye- Waller factor, however, this Debye- Waller factor has values much larger than would be 
calculated using simple models. A classical limit for incoherent scattering is obtained for high 
energies and temperatures. Between these limiting classical and quantum cases is another regime 
in which diffraction features appear that are broadened by the motion in the fast direction of the 
scattered beam but whose intensity is not governed by a Debye- Waller factor. All of these limits 
appear to be accessible within the range of currently available experimental conditions. 
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I. INTRODUCTION 



Recently, two laboratories have demonstrated that it is possible to diffract beams of 
fast-moving atoms and molecules from ordered surfaces under conditions of grazing incident 
angles, a method that has been given the name Grazing Incidence Fast Atom Diffraction 
(GIFAD).-^ 1 ^ The range of energies goes from several hundred eV up to tens of keV for the 
projectiles used which include He, H atoms and H 2 molecules. The target surfaces were 
single-crystal alkali halides. These energies correspond to de Broglie wavelengths several 
orders of magnitude smaller than the lattice spacing of the target crystal. Because of the 
long travel path near the surface one might expect sufficient dephasing to occur which 
would destroy all quantum coherence, and furthermore a simple application of the Debye- 
Waller attenuation due to thermal roughness of the surface would imply totally negligible 
diffraction intensities except possibly for the most grazing incidence conditions. However, 
when the incident beam was aligned along a high symmetry direction intense and clearly 
measurable diffraction peak signals were observed in the direction normal to the scattering 
plane. 

Two different theoretical groups have carried out numerical simulations of the diffraction 
process using propagation of wave packet techniques and have shown that one can obtain a 
reasonable description of the scattering intensities as functions of the initial beam param- 
eters such as energy and incident angle.— ^ These descriptions indicate that along the fast 
velocity component of the incident beam the interaction potential of the surface is averaged 
so that the total potential appears as if it were relatively smooth with one- dimensional cor- 
rugations in the high symmetry direction. Also, it has been shown that simple calculations 
of diffraction intensities based on the eikonal approximation for a one-dimensional repulsive 
hard wall potential can predict the diffraction intensities reasonably well. 1 

In this paper a theory is developed starting from semiclassical quantum mechanics that 
produces analytical expressions which describe the scattering process at levels of coherence 
ranging from purely quantum mechanical diffraction to the fully incoherent classical limit. 
The basic starting approximations are the eikonal approximation and assumption of a tar- 
get crystal whose thermal vibrations are in the harmonic limit. The scattering process is 
described as a series of sequential collisions along the path of the projectile as it moves near 
the crystal. Under these conditions very little energy is lost due to motion in the fast direc- 
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tion parallel to the surface and it is shown that the Debye- Waller factor becomes sufficiently 
large to permit measurable diffraction intensities. Under conditions where quantum peaks 
are observed the very small energy loss associated with the fast direction ensures that the 
dephasing is sufficiently weak to permit quantum scattering in the other two mutually per- 
pendicular directions. In addition to the purely quantum diffraction limit the theory can be 
taken to the classical limit through use of the Bohr correspondence principle. In between the 
clearly distinct quantum and classical limits there is a regime where the dephasing allows for 
distinct quantum peak features that are broadened by a limited coherence length and these 
quantum diffraction features do not obey Debye- Waller attenuation. Each of these regimes 
has its own signature characteristics that should be measurable in the energy, angular and 
temperature dependence of the measurements. 

In the remainder of this paper the theory is developed together with a description of 
the various possible quantum mechanical and classical limits in the following Section [Til 
In Sec. [TTTJ some numerically calculated results are shown and discussed. Conclusions are 
discussed in Sec. [TV] Some of the results used in developing the various models of Sec. [II] 
are developed in the Appendices. 

II. THEORY 

An appropriate starting point for developing the theory of scattering of atomic projectiles 
from a surface is the quantum mechanical transition rate for a particle of incident wave vector 
kj making a transition to the final state of wave vector ky given by 



where Tfi is the transition matrix element taken with respect to the final and initial states 
of the system of projectile plus target, Si and Sf are the initial and final global energies 
of the entire system. The angular brackets indicate an average over all initial states of the 
surface target and the Yl{n f } indicates a sum over all final states of the target. 

The usual procedure is to write the energy delta-function as its time-integral Fourier 
representation and to shift to the interaction picture in which the time dependence is gov- 
erned by the target hamiltonian^, a process that is sometimes called the Glauber- Van Hove 




(1) 
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transformation^: 



w(k f ,ki) = -2 
n 



+00 



dt e m - E f )t/h (T if (t)T fi {0)) 



(2) 



where Ef and E^ are the final and initial translational energies of the projectile atom. 

The fast atom diffraction experiment is one in which one component of the incident 
momentum, chosen as hki X , is much larger than the component in the direction normal to the 
surface which is chosen to be the z-direction. Thus, it will be a reasonable approximation to 
assume that the transition operator separates into the product of an operator in the normal 
direction and an operator in the directions parallel to the surface 
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The matrix element of the transition operator, taken with respect to final and initial states 
of the projectile, then separates into the product 
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*/(s,y) T xy ^i(x,y) 
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The squared normal transition operator 



<&i(z)) is proportional to the prob- 



ability of making a transition from ki Z to kf z . There are two useful limits in which this can 
be evaluated exactly. The first is the nearly trivial case of a surface whose potential is a rigid 
hard repulsive wall in which case it is a constant. This is the evaluation used in the stan- 
dard application of the eikonal approximation for elastic scattering.-^ The other useful limit 
is given by the classical motion of a grazing-angle projectile moving in a one-dimensional 
exponentially repulsive potential such as 
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If the atoms of the surface are subject to small random displacements about their equilibrium 
positions the scattered distribution becomes a lognormal distribution in the final angle^ 1 ^ 
as shown below in Appendix IA 21 
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where a 2 is the mean square displacement. The origin of the surface displacement is thermal 
motion, although the displacement appears static to the fast-moving projectiles because the 



collision times are much shorter than typical phonon periods. Thus a 2 is given by the mean 
square thermal displacement normal to the surface, which in the high temperature limit for 
a Debye model of the phonon density of states is 

where ks is Boltzmann's constant, T$ is the surface temperature, Mc is the surface mass 
and 0£> is the surface Debye temperature. 

The matrix elements of the operator T xy will be evaluated semiclassically within the 
eikonal approximation following the work of Bortolani and Levi^. For purely elastic scat- 
tering from a rigid potential the wave function in the asymptotic region far from the surface 
is of the form 

^( r ) = e *i* _ A(K) e i(Kl+K) - R e ik f* , (8) 

K 

where Kj is the surface-parallel component of the incident wave vector kj and K = Kj — Kj. 
The perpendicular component k/ z of k/ is determined by energy conservation Ef = E{ which 
implies 



"fz — ^iz 



K 2 - 2K • Ki . (9) 



The coefficient -A(K) of the outgoing asymptotic wave is related to the transition matrix 
element through 12 

(*/(*, y) %(x,y)) = ie iS f ^ A(K) , (10) 

where L is the quantization length. 

The eikonal approximation consists in first applying the Rayleigh ansatz, i.e., assuming 
that the wave function of Eq. (jSJ) is valid near the surface and then applying the correct 
hard wall boundary condition 

^(R,* = £(R)) = 0, (11) 

where £(R) is the corrugation function of the repulsive wall. The eikonal approximation is 
obtained by further assuming that kf z is weakly dependent on K which allows for a simple 
evaluation of the coefficient of the outgoing wave 

A(K) = f dR e- iK R e - iAk ^ R) , (12) 



where Ak z = k{ Z + kf z . 

Eq. (j!2p is the eikonal solution to the purely elastic problem. The more general problem of 
a potential with thermal vibrations can be treated by introducing the displacement function 
u(R, t) of the surface into the coefficient A(K) using the transformations R — ► R — u(R, t) 
and £(R) — > £(R)— w 2 (R, t). Applying this transformation to the transition matrix through 
Eq. ()12p casts the transition rate of Eq. (|2j) into the form 

w(k f M) = ± (^) 2 | (*/(*) \t\ *,(*)) f j^dte^-W (13) 

x _L y rfR y rfR' e -*K-(R-R') e -tf*,[€(R)-{CR')] ^ e iAk.u(R,t) e -iAk.u(R',0)^ 

where Ak = kf — k, = (K, Ak z ). Eq. fflBl is equivalent to the similar expression for surface 
scattering in the eikonal approximation developed in Ref. p^]. 

The thermal average in Eq. fflBl can be readily carried out in the harmonic approximation, 
and after taking proper account of the commutation relations of the displacement operators 
at different times and positions the result is 13 



w(k f ,k t ) = i, \(* f (z)\f z \Mz))\ 2 J ^ dte^-W (14) 

x J_ y dK J dR' e - iK -( R - R ') e -«*.K(R)-«M] e -2W(Ak) e 2W(Ak;R,R',t) _ 

The argument of the Debye- Waller factor exp{— 2W} is given by 

2W(Ak) = ((Ak-u) 2 ) , (15) 

and a standard approximation is to assume that all cross terms in Eq. (fl~5|) average to 
zero and that the mean square displacements in all three mutually perpendicular cartesian 
directions are the same. This leads to 

2 Wm ^^W =g£g. (16) 

The position and time dependent correlation function appearing in the exponential of 
Eq. ( JT41 is given by 

2W(Ak;R,R',t) = ( Ak ■ u(R, t) Ak • u(R, 0) ) . (17) 
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If the surface is periodic the correlation function also exhibits this periodicity and can be 
expanded in terms of the normal modes of vibration and appears as^ 1 ^ 

2W(Ak;R,R',*) = £ £ - » (18) 

Q,a'=l Q,u 

x e Q (Q,z/)e a ,(Q,z/)e jQ -( R - R ') {[n„(Q) + l] + n„(Q) e*^} , 

where e a (Q,z/) is the a cartesian component of the polarization vector of the mode with 
parallel wave vector Q and perpendicular quantum number v [y is discrete for surface modes 
and is a continuous variable for the bulk modes corresponding to the perpendicular wave 
vector) and N c is the number of modes. The mode frequency is u u (Q) and n u (Q) is the 
Bose-Einstein function 

MQ) = — • (19) 

e k B T s — 1 

Comparison of the correlation function of Eq. (fTTj) or (fT8]) and the Debye- Waller exponent 
of Eq. (fl5l) shows that 

W(Ak) = W(Ak;R = R',t = 0) , (20) 

a fact that becomes important later when discussing the classical multiphonon limit of the 
transition rate. 

Eq. ( p3| is the starting point for developing and discussing the various regimes of scat- 
tering that can occur in fast atom diffraction from surfaces. These regimes range from the 
purely quantum mechanical case to completely incoherent classical scattering and in the 
subsections below we discuss several of these possibilities. 

A. Classical Scattering 

The classical regime occurs when the scattering is completely incoherent as for example 
when a large number of phonons are transferred in the collision and the quantum coherence 
of the incident beam is completely destroyed. One case in which this occurs is when the 
Debye- Waller factor becomes negligibly small implying that its argument 2W is large.— »^ 
The value of the Debye- Waller agrument is essentially a measure of the number of phonons 
transferred in the collision, and it will be large when either Ak or the mean square displace- 
ment (or effectively the temperature T5) becomes large. When the Debye- Waller argument 
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2W becomes large the correlation function 2W of Eq. (ITS]) also becomes large and further- 
more its most important contributions come from the region of small t and small R — R'. 
In this case the leading terms in the expansion of 2W are 

2W(Ak; R, R', t) « 2W(Ak) (21) 

3 h 

- it Ak *&K 2N M e a (Q,u) e a ,(Q,v) 

a,a'=l Q.f ^ ^ 

- t 2 J2 Ak " Ak '* E ^^e Q (Q,z,)e a ,(Q,z.) [2n,(Q) + 1] 

a,a'=l Q,u C C 

- £ Ak * Ak « E h[ WMuJQ) e °(Q'") e °'(Q,") [2n,(Q) + l] , 

where the Debye- Waller argument 2W(Ak) is the same as that of Eq. (fT5l) 

The simplest approximation for evaluating the expansion of Eq. (1211) is to assume a 
Debye-Model with bulk-like symmetry in which all cross terms vanish and for this case it 
simplifies to 



2W(Ak; R, R, t) « 2W^(Ak) - iAE t - AE k b T s t 2 S \ „ '- , (22) 

2h v R 

where the classical recoil energy is given by AE = h 2 Ak 2 /2Mc, and the velocity parameter 
vr is given by the relation^ 

T = E A *« A fe E o^ Q '^ 6a(Q '^ e "'(Q^) [2^(Q) + l] , (23) 

where R is a unit vector in the direction of R — R'. The parameter vr of Eq. (1231) is a 
weighted average of phonon speeds parallel to the surface and it can be evaluated for several 
simple models of the surface phonon density^. These models produce values that are of 
order of the bulk acoustic phonon velocities or the Rayleigh wave velocity. 

With the simple evaluation of the correlation function of Eq. ( |22l) the transition rate 
of Eq. (|T3l) can be evaluated, a process equivalent to a steepest descents approximation. 
Ignoring for the moment the final term involving the spatial dependence in Eq. ( 122|) the 
result is 

»(k,,kO = f (^) 2 |(*,(z)|t,| *,(*)) f \A(K)\ 2 (24) 
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Eq. (1!M1) is an expression that describes a single classical collision. For example, if the 
prefactor transition matrix elements are chosen to be a constant, i.e., if 



the choice appropriate for hard sphere scattering where the classical cross section is uniform 
in all angular directions, then it becomes the classical transition rate well-known in neutron 
scattering^. It is also an expression that has been demonstrated to explain single collision 
events in low energy alkali ion scattering from metal surfaces^ 1 ^ as well as energy-resolved 
measurements of rare gas scattering from molten metal surfaces. 20 

The spatial dependence appearing in Eq. (1221) that was ignored in deriving the transition 
rate of Eq. (1241) will become of importance to the discussion below in Sec. Ill Dl However, it is 
of interest to note here that another classical expression for the transition rate was obtained 
by Brako and Newns under the assumption of a smooth surface with thermal vibrations. 
If all terms of Eq. ( 1221) . including the term in the spatial dependence, are used to evaluate 
Eq. (1T31) . but at the same time assuming that the surface is flat, i.e., setting the corrugation 
function £(R) equal to a constant, the result is^ 



where the transition matrix prefactors have again been set equal to a constant. This expres- 
sion has also been shown to describe rare gas scattering from metal surfaces.— ^ For the 
conditions considered here for fast atom scattering at grazing incidence conditions Eq. fT26|) 
will not be a useful approximation. For the expected values of v r, which is in the neighbor- 
hood of phonon speeds, the Gaussian-like function in parallel momentum will be very broad 
and will not have an appreciable effect on the overall transition rate. 

Thus, the basis for describing classical scattering will be Eq. (124")) . but as mentioned 
above it describes a single scattering event in which the projectile particle is assumed to 
make a transition from its initial state of momentum frkj to the final state of momentum 
hhf. However, the discussion in Appendix [Bl shows that this is not at all the case. Instead a 
grazing angle projectile will follow a trajectory along the surface and will encounter succes- 
sive, incoherent collisions with approximately N atoms before being deflected into the final 
state of momentum Kkf. 




(25) 
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Under grazing angle conditions, which imply small total scattering angles 9, the transition 
rate of Eq. (l2lj) becomes a very sharp Gaussian as a function of final energy Ef with an 
average energy loss given by 

Ei - E f « AE « /i^0 2 , (27) 

where /i is the mass ratio and the final term on the right of Eq. (127]) arises because in a 
small angle collision with little energy loss |Ak| This same result can be obtained 

from the Baule conditions for scattering through small angles as shown in Appendix [B] 

As shown in more detail in Appendix [B] the multiple scattering events along the trajec- 
tory can be viewed as a series of N convolutions of the Gaussian-like differential reflection 
function. The energy loss during these N collisions is N times that of Eq. (127)) but with the 
total scattering angle 9 replaced by the fractional angle suffered at each of the N collisions 
which is approximately 9/N. Given that the square of the full width at half maximum 
(FWHM) of the Gaussian approximation to Eq. (12^)) is proportional to AksTsAEo and that 
for the convolution of Gaussians the squared width of the convolution is the sum of the 
squares of the individual widths, the final result is a transition rate that looks exactly like 
Eq. fl2M except that the recoil energy is replaced by 

AE 

AE ~ ~JT ~ ~W~ ■ (28) 

As shown in Appendix [B] this energy loss is very small and, moreover, is nearly entirely 
associated with the motion in the slow directions. The energy loss in the fast direction is 
negligible in comparison because it can be readily shown that Ak 2 x is smaller than Ak\ + Ak 2 
by a factor of 9} which for grazing angles gives a reduction of several orders of magnitude. 
This small energy loss in the fast direction implies that the scattered projectiles will lie close 
to the energy conservation circle defined by 

kj z + kf y = k iz + k iy (29) 

which also implies that the total scattering angle 9 is given in all cases to a good approxi- 
mation by 29 i. 

Combining the above simple approximation of the total scattering angle with the eval- 
uation of N = 6/Ta9i for the exponentially repulsive potential of Eq. (j3J) as obtained in 
Appendix [B] leads to 

AE = = 2 m 
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where a is the distance between surface atoms in the fast direction. 

This then leads to a final expression for the classical transition rate which looks similar 
to Eq. f)24p but with important modifications: (1) AE is replaced by the much smaller AE, 
(2) the scattering amplitude |v4(K)| 2 separates into the product of |v4(Afca.)| 2 \A(Ak y )\ 2 . 
However, in the fast direction the surface is smoothly averaged and |A(A/c x )| 2 is a constant. 
For the slow direction parallel to the surface ^(AA;,,)! 2 is the probability of being deflected 
out of the scattering plane. In the classical limit it will contain the rainbow features of the 
surface profile, which for example, can be calculated from the classical limit of the eikonal 
approximation of Eq. (fT2l . (3) Finally, the transition matrix element in the normal direction 
is given by the lognormal distribution of Eq. (1A22P with Of given by the relation for the 
energy conservation circle of Eq. ( |29i) . 

The relationship between the transition rate and the differential reflection coefficient actually 
measured in the experiments is made explicit in Appendix O 

The differential reflection coefficient described by Eq. (13TT) has a width in the final po- 
lar angle (or equivalently in the direction of kf z ) dictated by the width of the lognormal 
distribution P(9f). Its spread in the azimuthal direction (or equivalently the kf y direction) 
is determined by the form factor |y4(Afc y )| 2 which depends on the shape of the surface po- 
tential along the y-direction. The intensity will exhibit peaks at the rainbow angles which 
correspond to the inflection points in the surface profile and the spread in the azimuthal 
direction will be limited to the classically allowed region, which is typically within the limits 
of the rainbow features allowed by the largest tilt angle of the surface profile. 



B. Purely Quantum Mechanical Scattering 



The case of purely quantum mechanical scattering in all three dimensions is formally no 
different from the case of ordinary atom-surface diffraction at low energies. This is most eas- 
ily discussed in terms of the eikonal approximation 11 or the kinematical approximation^ in 



which case the transition rate becomes the same as Eq. (IT4I) except with 
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taken to be constant. 

^-h&i !Z ***-*™ (32) 

X _L J dK J dR' e - lK -( R - R ') e -^K(R.)-?(R-')] e -2^(Ak) e 2W(Ak;R,R',i) _ 

This is then usually developed into a series ordered in numbers of phonons transferred by 
expanding the exponential of the correlation function. The elastic scattering is the zeroth 
order term in this series given by 



X 



' f dR e~ iKR e"^( R ) 2 



L 2 



e -2W(Ak) 



The integral over time becomes a simple S(Ef-Ei) stating the conservation of energy. If the 
corrugation function £(R) is periodic in both directions parallel to the surface the eikonal 
approximation matrix element reduces to an integral over a single unit cell multiplied by a 
Kronecker 5-function restricting the parallel momentum transfer to reciprocal lattice vectors 
G of the surface 

w(k f X) = y (^tVE \A(G)f e-™^ 8 k , g 5(E f - Ei ) , (34) 

V m J Q 

where the scattering amplitude is 

A (G) = -L / dR e~ lG R e" iAfc ^ R ) . (35) 

Juc 

As stated above, Eq. (|34j) is the eikonal approximation result for diffraction from a two- 
dimensional periodic and strongly repulsive surface. However, in general if h 2 kf z A(G)/mL 
is replaced by the transition matrix element of the full elastic interaction potential Eq. (|34|) 
becomes exact and is not an approximation. 

There are several reasons why it is unlikely that the purely quantum mechanical diffrac- 
tion described by the case of Eq. (1341) . other than possibly a peak at the specular position, 
will be observed in the kf x direction in fast atom diffraction: (1) for energies of hundreds of 
eV the Debye- Waller factor would be be extremely small except for the most grazing inci- 
dent conditions, while experimentally it is known that diffraction features in the transverse 
fc/ y -direction are observed for energies into the keV range at incident angles of the order of 
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one degree above the surface plane.- 1 ^ (2) The experimental evidence indicates that in the 
fast direction the atomic projectile experiences an averaged potential that is not corrugated, 
implying that the matrix elements will be small for non-zero order diffraction peaks in the 
fc^-direction.— (3) The conditions of energy and parallel momentum conservation determine 
the value of the perpendicular momentum and similar to Eq. (jH]) above it is restricted to the 
Laue circle given by 

k% = kj z - G 2 - 2G Kj . (36) 

With ki X » k iz the cross term on the right hand side of of Eq. (1361) will dominate and 
for positive G x values kj z < and those beams will be evanescent.- For negative G x values 
kf z will become substantially larger than ki Z implying that the final diffraction angle will be 
larger than 9i and will not satisfy grazing conditions for the exiting beam. This again leads 
to the conclusion that such diffraction features with G x ^ would have negligible intensity. 

C. Quantum Behavior in only Two Dimensions 

A more likely scenario for producing diffraction is that the problem reduces to a mixed 
quantum-classical case, with purely quantum mechanical behavior in the two slow directions 
(the z- and ^/-directions) and classical behavior in the fast direction. In this case the motion 
of the atomic projectile in the fast direction is regarded as a series of N incoherent and 
successive collisions with the surface atoms just as in Section III Al above meaning that 
the lognormal distribution in final angles is retained. However, the quantum treatment of 
Section III Bl above is carried out for the y-direction. The result is 

-(k/,k 4 ) = | (^f Y P(9 f ) \A(G y )f e-^' N 5 Aky , Gy P (E f ) , (37) 

Gy 

where p{Ef) is the density of final states in translational energy identified in Sec. Ill Dl below. 
The scattering amplitude in the y-direction calculated with the eikonal approximation is now 

A(G y ) = — [ " dy e~ iG y y e" iAfc ^ (s/) , (38) 

a y JO 

which for a simple one- dimensional sinusoidal corrugation function £(y) = ha y sin (^yj 
where h is the corrugation strength leads to the well-known Bessel function solution 8 

A(G y ) = J lm] (ha y Ak z ) , (39) 
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where J n (q) is the ordinary Bessel function of order n and argument q. 

An important modification has been applied to the Debye- Waller factor, it appears di- 
vided by a factor of N. If the interaction with the surface is considered to be a series of 
incoherent collisions with N surface sites, a Debye- Waller factor should be associated with 
each of these collisions. Such a multiplicative cumulation of Debye- Waller factors is used 
in LEED for describing the multiple collisions that the electrons incur as they penetrate 
the surface.— However, the Debye- Waller argument 2W for each of these small collisions 
is smaller than that of Eq. (|T6j) by a factor of 1 / iV 2 for reasons similar to those that were 
discussed in connection with Eq. ( 1281) above, i.e., because |Ak| m ki9 and for each small 
collision 9 m 28i/N. Thus, this argument implies that the Debye- Waller factor of Eq. (l37p 
is transformed according to 

e -2W _^ (e-2W/N^ = e -2W/N (4Q) 

This same result can be arrived at by a somewhat different argument .— The Debye- Waller 
factor is the average of the phase induced by the thermal displacement 

e -W = ^Ak-u) = e -|((AM 2 ) . (41) 

If the total phase is considered to be the sum of the individual phases picked up in each of 
the small collisions 

N 

Ak-u ^ ^ Ak n ■ u n , (42) 

n=l 

and if the individual scattering events are considered uncorrelated the sum of all cross terms 
in the mean square thermal average ((Ak ■ u) 2 ) will go to zero leading to a sum containing 
N mean square phases. Assuming a similar phase for each small collision again results in 
the same expression as in Eq. (j40]) . 

There is a seemingly apparent conflict inherent in the result of Eq. (1371) because the 
density of states in final energy p(Ef), if evaluated strictly, is a 5-function in the difference 
between final and initial translational energy. However, if the motion in the fast direction 
is classical then the translational energy of the projectile cannot be conserved because some 
energy will be exchanged with the surface through the classical collisions. The conservation 
of energy is very nearly preserved because the exchange of energy associated with the fast 
direction is negligibly small for the same reasons as explained above in Section [II Al i.e., 
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as long as all angles are grazing Akf. is smaller than the momentum transfer in the other 
two directions by a factor of order Of. Thus the condition of energy conservation is weakly 
relaxed, and this allows for the range of kf z values implied by final polar angle width of the 
lognormal distribution. The expression for p{Ef) that allows for a small exchange of energy 
is determined in Eq. (1521) of Section III Dl below. 

The scattered intensity spectrum described by Eq. (|37|) exhibits a number of interesting 
features. In the kf y direction there are narrow diffraction peaks whose intensity are given by 
the transition matrix element of the potential which in the eikonal approximation is Eq. (1381) . 
In the actual experiment, the widths in this direction would be limited by the transverse 
coherence length of the experiment. In the kf z direction these diffraction features are broad 
with a width determined by the angular spread of the lognormal distribution. All of these 
diffraction features lie roughly along the Laue circle given by 

k% = tfz ~ Gy — 2G y kiy , (43) 

or equivalently 

cos# f = cosflj • (44) 

1 COS0/ v ; 

In each of the diffraction features kf z increases from the bottom to the top of the broad 
streak which means that the intensity can vary considerably as is seen from the phase 
variation in the scattering amplitude of Eq. (|38|) . This quantum interference effect can be 
either destructive or constructive, and when it is destructive it can even make portions of 
the streak disappear. This quantum interference effect can also cause a given streak to 
appear as if it were shifted above or below the Laue circle. Finally, the Debye- Waller factor 
is considerably larger than that implied by the simple expression of Eq. (1T61) making it 
possible to observe diffraction peaks with measurable intensities. The Debye- Waller factor 
essentially depends on momentum transfer only in the quantum mechanical directions, i.e., 
Ak 2 w Aky + A/c 2 because at such small scattering angles A/c 2 is negligible in comparison. 
All of these features are observed in the experimental measurements^. 



D. Near- Classical Scattering with Spatial Correlations 



The two cases above in Sees. Ill Bl and lll Cl discuss conditions in which quantum mechanical 
diffraction peaks can arise in the scattered intensity. There is yet a third case in which spatial 
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i /fr 2 L \ 2 r +oc 

;(k/ ' kj) = ¥ V^l) P{6f) J dtei{Ei ' EfW% ( 45 ) 



correlations will give rise to diffraction features in an otherwise nearly classical scattered 
spectrum. This case is similar to the case of diffraction features observed in the multiphonon 
background in electron diffraction experiments such as RHEED where they are known as 
Kikuchi lines.— 

The starting point for this treatment is Eq. (j!4p assuming an averaged flat potential in 
the fast direction and a corrugation function £(y) of periodic length a y in the ^-direction 

w{ 

L 2 

The periodicity of the corrugation function is expressed as £(y+na y ) = £(?/) with n an integer 
and the fact that the correlation function W depends only on the difference in displacement 
y — y' allows one to replace the spatial integrals by a sum of integrals over a single unit cell. 
Since the correlation function has very weak spatial dependence over the distance of a single 
period (or equivalently, under the assumption that each unit cell vibrates rigidly) this leads 
to 

i /h 2 h \ 2 r + °° 
Mk /)kj ) = P{fi f ) J ^ dte«*-*fM (46) 
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Making the classical limit expansion of the correlation function as in Eqs. (T2TI) and (|22|) 
leads to 

-.2 1. \ 2 



W 



[kfM) = ^ (^) P ( d f) dt e^- E ^ h e- tAEt ~ AEk * Tst2 (47) 

+oo +oo 

X \A(Ak y )\ 2 ^ e~ iAkyavin ' n ' } e -^ Ek bT S {n-n')^al/2h 2 v% _ 



n=— oo n'=— oo 

where Ai?o has been replaced by A.E of Eq. fl30l) to account for the multiple collisions with 
iV surface unit cells along the fast direction, and as in Eq. (|38|) 



1 



A(Ak y ) = — / e" iAfc ^ e - iAfc ^ } . (48) 
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The discrete sums reduce to a single summation given by 



dy 

w n=— oo 



^ e -i&k y a y n e -AEk b T s n 2 a y /(2h 2 v 2 R ) ^ 



A ^» e 2h 2 v 2 R (Ak y -G y ) 2 /(4k B T s AE) 

a y y/2irAEk B T s ^ 



where L is a quantization length parallel to the surface. Clearly, this summation is a periodic 
function in Ak y — ► Ak y + G y where G y = 2'Kv/a y with v an integer, and the right hand 
side of Eq. ( 119]) becomes a good approximation if AEkBTsa y / (2h 2 vj t ) < 1. 

The time integral in Eq. f|4"5]) is the same Gaussian integral encountered in Section III Al 
above so the final result for the transition rate is 

x MfAA; )l 2 V e 2^4( Afc !/-G ! /) 2 /(4fc fl T s Ai?) 

1 1 ^ V2nAEk B T s 

Eq. (|50|) describes periodic diffraction features lying roughly centered on the Laue circle 
of Eq. (143p . In the direction these features are Gaussian- like with a width parameter ko 
given by 



<&Ek b T s al 

-2iNT ■ (51) 

In the kf z direction they appear as long streaks whose width is governed by the lognormal 
distribution in final polar angles. The strength of these peaks at any point (kf z , kf y ) in the 
wave vector space perpendicular to the fast direction is given by the form factor |A( Afc y ) | 2 
which in this treatment is approximated by the eikonal expression of Eq. (|48|) . This form 
factor can vary substantially over the kf z extent of the diffraction feature, just as discussed 
above in connection with Section III CI even to the point of making parts of a diffraction 
feature disappear because of destructive quantum interference arising from the corrugation 
within the unit cell of the potential. There is no Debye- Waller factor associated with these 
diffraction features. Instead the temperature dependence varies inversely with a power of 

—3/2 

the temperature, in this case T s , typical of classical scattering. Also typical of classical 
scattering, the widths of all peaks increases as y/Ts for temperatures large enough where 
zero point motion is not important. The energy spread within these peaks is quite narrow 
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and given by 

(E ) 1 ex ( (Ef-Ej + AE) 2 

P[ f) ' sJ^AEk B T s I AnAEk B T s 

which approaches the energy ^-function S (Ef — Ei + AE) for small AE. Eq. (|52|) also 
identifies the density of final energy states p (Ef) introduced above in Eq. fl37|) . 




III. CALCULATIONS 



In the following a few calculations are presented to illustrate some of the different cases 
presented in the sections above. The system considered is 3 He scattering from LiF(OOl) 
along the (110) direction with an incident energy of 200 eV and an incident polar angle of 2° 
with respect to the surface. With the LiF lattice vector of b = 4.03 A the effective reciprocal 
lattice vector is G y = 2n\/2/b = 2.21 A -1 . 

For the calculations the surface Debye temperature is chosen to be 530 K (a value cor- 
responding to a bulk Debye temperature of 750 K divided by y/2), the mass M c = 19 
amu corresponding to that of a fluorine atom, the exponential potential range parameter is 
T = 2.0 A -1 , the surface temperature is 300 K and the velocity parameter is chosen to be 
vr = 4000 m/s which is the Rayleigh wave speed of LiF(OOl).— The scattering amplitude 
A(Ak y ) that determines the relative intensities of the diffraction features is chosen to be the 
eikonal Bessel function approximation of Eq. (139]) for a simple one-dimensional sinusoidal 
corrugation function with the arbitrary choice h = 0.023. The evaluation of the number 
iV = 6/TaOi of LiF lattice sites encountered along the path of the He atom during the 
interaction is approximately iV ~ 21 for these conditions. 

Fig. [1] shows a graph plotting the differential reflection coefficient dR/dOfdifif as a function 
of 4>f (measured in degrees) with the polar angle chosen to lie along the Laue circle of 
Eq. (1431) or (jUJ). The intensity plotted on the vertical axis is measured in arbitrary units. 

The broad diffraction features centered about the diffraction points are calculations of 
the nearly classical transition rate of Eq. (1501 . The sharp vertical lines located precisely 
at the diffraction positions Ak y = G y are the quantum diffraction peaks calculated from 
Eq. fl37|) . The relative intensity of the two sets of calculations with respect to each other is 
arbitrarily chosen. In this particular example the zeroth order or specular diffraction peak 
at cf)f = has nearly zero intensity and peaks of up to order 4 are visible on either side of 
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the specular position. 

Fig. [2] shows the diffraction feature calculated with the nearly classical transition rate of 
Eq. (1501) but this time plotted as a function of 9f measured in degrees for values of 0/ fixed 
at the positions of the first and third order diffraction peaks in the upper and lower panels, 
respectively. This shows a rather broad feature with a FWHM of about half a degree which 
is essentially the width of the lognormal distribution in final polar angle. The polar angle 
distribution for the diffraction peaks calculated with the quantum model of Eq. fl37|) would 
be identical. Even for the third order diffraction peak a small shift of the distribution along 
the Laue circle to polar angles smaller than Of = 2° is evident. 

Figs. H] and [2] show the intensity integrated over all final energies because this is what 
is experimentally measured by the channel plate detector. In order to illustrate the energy 
dependence the full three-dimensional differential reflection coefficient dR(hf,]<.i)/dEfdQf 
is plotted as a function of final energy Ef in Fig. [3] for final angles at the specular position 
and at the second order diffraction peak in the upper and lower panels, respectively. The 
calculations are carried out with the nearly classical transition rate of Eq. (1501) . The two 
peaks are quite similar, and in fact the energy spectrum in the vicinity of all diffraction 
peaks looks very much the same. This figure demonstrates the remarkably small energy 
losses experienced by these projectiles in spite of the high incident energy and relatively 
large incident angle. The full width at half maximum is less than 50 meV and the most 
probable final energy is approximately 10 meV smaller than the 200 eV incident energy, or 
a fractional energy loss of about 5 x 10~ 5 . This energy loss is in agreement with the value 
of about 7.2 meV given by the expression of Eq. f l30l ). 

IV. CONCLUSIONS 

The objective of this paper is to discuss the variety of scattering regimes that can occur in 
the diffraction of fast atomic and molecular projectiles at grazing angles from surfaces. Such 
a discussion should begin with mention of purely quantum diffraction in all three spatial 
dimensions. This would be essentially no different than a description of the well-known 
technique of ordinary He atom diffraction at low (thermal) energies, and is unlikely to apply 
to the present experiments because of the high energies and correspondingly extremely large 
velocities in the fast direction. 
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However, here we review the salient features of three regimes of diffraction and scattering 
that appear to be possible. These are (A) quantum behavior in two dimensions, (B) near- 
classical scattering with spatial correlations, and in addition, (C) possible contributions to 
the incoherent background are briefly reviewed. Each of these regimes has very characteristic 
features in their dependence on the experimentally controllable variables. 

A. Quantum behavior in two dimensions 

In this case, discussed in detail above in Sec. Ill C[ the motion in the fast direction, parallel 
to the surface, is classical while the motion normal to the surface and also that parallel to the 
surface but transverse to the fast direction is quantum mechanical. In spite of the classical 
motion with multiquantum energy transfer in one of the three dimensions, quantum features 
can be observed in the other dimensions because the energy losses in the fast direction are 
exceedingly small, too small to cause a complete loss of quantum coherence. 

This regime is characterized by diffraction features that appear as streaks in the /in- 
direction, with the length of the streaks governed by the lognormal distribution P{9f) in 
final polar angles Of. In the /c/^-direction these peaks are narrow, their widths being limited 
by the instrument resolution. These peaks are observed to be distributed roughly around 
the Laue circle whose wave vector components are given by Eq. (I5ET) . 

The intensity dependence of these diffraction features obeys a variant of the classic Debye- 
Waller behavior. The temperature dependence of the peak intensity is an exponentially 
decreasing function of temperature divided by a power of Tg, and in the treatment shown in 
Eq. (1371) it is seen that power is one. The argument of the Debye- Waller factor is considerably 
reduced over the standard value that would be calculated according to Eq. ffT6]) using the 
initial and final momenta of the diffraction features. The 2W argument of the Debye- 
Waller factor is effectively divided by N, the number of collision sites along the path of 
motion. It is this effect that allows the Debye- Waller factor to be large enough that the 
diffraction features are measurable. Also, the motion in the fast direction does not contribute 
appreciably to the momentum transfer in the Debye- Waller argument of Eq. (Tl6l) and this 
is the cause of the extremely small energy losses in that dimension. For beams incident 
along high symmetry directions of the surface, where the Laue circle of Eq. ( 156]) appears 
on the channel plate detector as symmetric about the specular spot, it is of interest to 
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note that the Debye- Waller exponent 2W gets smaller, and hence the Debye- Waller factor 
becomes larger, with increasing G y reciprocal lattice vectors, and this is because of the 
associated decrease of the final perpendicular momentum transfer. If the incident beam 
is slightly miss-aligned with respect to a high symmetry direction, which would imply a 
small non- vanishing Kj in Eq. (I36p . the Laue circle becomes skewed with respect to the 
specular point and the preceeding statement becomes true only for those reciprocal lattice 
vectors G y whose corresponding perpendicular momentum transfers decrease, since some 
perpendicular momentum transfers may increase. The temperature dependence of the width 
of the diffraction features differs in one respect from the standard Debye- Waller behavior. 
In the /c/y-direction, the widths are independent of temperature, which is the normal D-W 
behavior. However, the widths of the streaks in the fc^-direction increases in proportion 
to \fTg. This increase in width comes from the lognormal distribution in Eq. (1371) which 
defines the overall length of the streaks. It is also of interest to note that the lognormal 
distribution shows that the width of the streaks in the fc/ 2 -direction will be proportional to 
T, the range parameter of the exponentially decreasing interaction potential in the surface- 
normal direction. 

The relative intensities of the diffraction features associated with different G y reciprocal 
lattice vectors, as well as the relative intensity within the streak of the diffraction feature due 
to a single given G y , is governed by the corrugation landscape of the potential as expressed 
in the transition matrix amplitudes A(G y ) of Eq. (!37l) . For example, the corrugation of the 
surface can produce rainbow effects and supernumerary rainbow effects in the diffraction 
peak distributions, as has been recently reported 3 . Intensity modulation due to the shape of 
the potential can also be observed within the streak of a single diffraction feature and such 
effects can cause the streak to appear "shifted" above or below its expected position on the 
Laue circle, as is apparent in the observations-. 

In addition to the above responses of the diffraction features to the controllable exper- 
imental variables, it is of interest to discuss several possibilities that might be open for 
future experimental investigation. The first of these is the unlikely possibility of observing 
quantum diffraction effects in the high velocity direction, due to reciprocal lattice vectors 
in the G^-direction. It is of interest to note the reasons why this is unlikely. For the case 
of positive G x (meaning the direction of G x is parallel to the fast velocity component) at 
the energies used in the present experiments the corresponding kf z components given by 
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the Laue equation of (1361) become rapidly evanescent and hence unobservable. On the other 
hand, for negative G x the value of kf z becomes sufficiently large that the diffraction fea- 
tures are projected several degrees in polar angle above the specular position. This means 
that these peaks would also be unobservable because a flat ordered surface cannot scatter 
intensity into such large scattering angles. 

A possible effect closely related to supernumerary rainbows is the observation of quantum 
interference due to scattering of terraces at different atomic height intervals on the surface as 
seen in thermal energy He atom diffraction^ and in reflection high energy electron diffrac- 
tion^. This phenomenon should also manifest itself in the present fast atom diffraction-type 
experiments, although the effect may be limited by the relatively short coherence length of 
the experiment. 

An effect observed very early in the history of He atom diffraction at thermal energies 
was the observation of selective adsorption, or resonances with the bound states of the 
potentiaP^. In principle these resonances, which are characterized by sharp increases 
or decreases in diffraction peak intensities as an experimentally controllable parameter is 
varied through resonance conditions, should be present in the current experiments. The 
simplest statement of the condition for selective adsorption resonance is that the value of 
h 2 k'j z /2m must be negative (corresponding to an evanescent beam) and nearly equal to the 
energy of a bound state of the potential. In the present experiments, the simplest way 
to look for selection adsorption events would be by slightly varying the incident azimuthal 
angle fa away from a high symmetry direction. When the incident beam is aligned along 
a high symmetry direction with fa = 0° the diffraction features observed on the channel 
plate detector lie roughly along the Laue circle of Eq. (136|) which appears symmetric with 
respect to the specular spot. However, if fa deviates slightly from zero, for example upon a 
small azimuthal rotation of the target crystal, the observed Laue circle becomes skewed and 
asymmetric. Along with this asymmetric shifting of the Laue circle, many of the low-index 
diffraction features will become evanescent and have kf z values that pass through conditions 
for resonance. 

One of the most important uses of thermal energy He atom diffraction has been the 
measurement of single phonon transfers and surface phonon dispersion relations. Because 
of the large disparities of energy between that of typical phonons and the kinetic energies 
used in theses experiments it is unlikely that most single surface phonon inelastic events will 
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be measurable. However, it is worth mentioning that there is one type of inelastic process 
that could be measured by these experiments, and this would be energy exchange with 
dispersionless excitations, e.g., surface phonons, surface plasmons, surface polaritons, etc. 
whose energies do not depend on momentum transfer parallel to the surface. For example, 
diffraction should still be observable even in the event that the atomic projectiles exchange 
one or more dispersionless inelastic quanta because such processes would not destroy the 
quantum coherence of the system. 

B. Near Classical Scattering with Spatial Correlations 

This case is described in detail in Sec. Ill Dl and is analogous to Kikuchi lines observed 
in RHEED due to simultaneous multiple phonon transfer and diffraction. The diffraction 
features appear similar to those due to two-dimensional diffraction described in the above 
section, except that there is no Debye- Waller factor and there is temperature dependent 
broadening in both the kf z and kf y directions. This type of diffraction could appear under 
conditions where the 2-D quantum behavior of the above section has negligible intensity, or 
the two can appear simultaneously as shown in the calculations presented in Fig. HJ The 
expression for the transition rate is given in Eq. (1501) and the major difference with the 2-D 
diffraction case of Eq. (j3"7|) is that the product of the Debye- Waller factor and the discrete 
delta-function dependence in parallel momentum transfer have been replaced by a periodic 
density of states consisting of broader peaks. The physical reason for these differences is 
the fact that multiple phonon quanta are transferred in this case, and the spatial region of 
coherence from which an interference diffraction pattern is built is no longer infinitely large, 
but becomes finite in extent. 

The diffraction features are still manifest as relatively long streaks in the fc^-direction 
whose length is governed by the lognormal distribution and these streaks appear aligned 
roughly along the Laue circle. In the /cjy-direction the peaks are now broadened beyond the 
instrumental limit, and the temperature dependence of the FWHM of this broadening is 
the same in both directions, i.e., proportional to v^Ts- The temperature dependence of the 
maximum intensity of the diffraction features, as shown in Eq. ( l50l) . is inversely proportional 
to T 3 / 2 . 
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C. Classical Incoherent Limit 



The classical incoherent limit is essentially a part of the background intensity that should 
be observed in the experiments. In this paper two different mechanisms from which it can 
arise are discussed. The first of these, presented in Sec. Ill A\ arises from scattering from 
regions of the surface whose heights are randomly distributed but whose average interac- 
tion potential is a repulsive exponential with range parameter T. The appearance of this 
background intensity is a feature broadened in the /^-direction by the lognormal distribu- 
tion, and broad in the /c/^-direction according to the extent permitted by the associated 
classical scattering form factor ^(AA^)! 2 . The temperature dependence of the maximum 
of this background intensity varies approximately inversely with the surface temperature, 
while the FWHM is proportional to y/Tg. Because the peak width in the polar angle Of is 
governed mainly by the lognormal function, this width is only weakly dependent on inci- 
dent energy. The peak intensity also is inversely proportional to T^/JI, while the FWHM is 
directly proportional to this same quantity. 

The total energy loss is evaluated in Appendix |X] in the classical limit of large incident 
energy and high temperatures and is given by AE = 2fiEiTa8f /3. This energy loss, for 
the present grazing-angle geometry with small 9i is negligible and roughly five orders of 
magnitude smaller than the incident energy. Most of this energy loss can be ascribed to 
motion in the direction perpendicular to the surface, and it can be readily shown that 
the the fraction of energy that can be ascribed to diminishing the fast hkf x component of 
momentum is smaller than the above AE by a factor proportional to Q\. It is of interest to 
note that AE is proportional to the range parameter of the repulsive exponential potential, 
implying that measurement of the energy losses over a range of incident energies, or when 
combined with measurements of intensity and FWHM dependence, may provide a means to 
measure this important physical parameter. 

The other, and radically different, type of background expected to be measured in high 
energy atomic scattering experiments is that due to binary collisions with topological defects 
such as adatoms or step edges on the surface. This contribution is described by Eqs. (|2~4"|) 
taken together with the constant form factor of (1251) . This is discussed further starting with 
Eq. (IB2I) of Appendix [Bl This type of binary scattering is spread over a broad range of 
energies and angles, with a FWHM proportional to \fE{Ts and an intensity at its peak 
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maximum inversely proportional to the same quantity, and its general features under the 
assumption of hard-sphere scattering are outlined in Eqs. (IB3I) to (IB 71) . In addition to 
its larger energy and angular widths, this type of background can be distinguished from 
that discussed in the two paragraphs above by its large energy losses, the fractional energy 
loss under the grazing angle conditions of the present experiments would be approximately 
fiEi8? as opposed to the cubic behavior in 8i discussed above for AE. This type of classical 
incoherent scattering is well-documented and has been extensively measured in experiments 
in low-energy ion scattering over energy ranges similar to those of interest here^. 
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APPENDIX A: CLASSICAL MOTION IN AN EXPONENTIAL POTENTIAL 



For a point particle moving in a one-dimensional exponentially repulsive potential 

V(z) = V e~ Tz , (Al) 
the expression of conservation of energy is 

E z = l -mv 2 z + V e~ Vz , (A2) 
where E z is the energy associated with the normal motion and the normal velocity is 



*<*> = § - V^ (1 -^ 0,) ' (A3) 

where zq = ]h(E z /Vq)/T is the classical turning point. Equation (1A3[) admits an analytic 



solution that can be parametrically expressed in terms of Vi Z = \j2E z /m as 

2 (\ + e~ Tvizt \ 
z(t) = z + v lz t + - lnf J , (A4) 

as can be readily verified by direct differentiation. 



1. Energy loss in a classical grazing incidence collision 

Beginning from the known expression for the energy transfer in a two-body collision with 
a small total scattering angle 9 

AE = Ei-Ef = ^9 2 Ei , (A5) 

where /i is the mass ratio, the small fractional energy loss occurring over a short distance 
parallel to the surface during which the particle is deflected by the small angle 56 is 

V 

The assumption of grazing angle scattering also implies that the parallel velocity v x » v z 
implying that 

56 5v z {t) 1 5 2 z{t) 1 



5E = nWEi = fiEA S 4) (5t) 2 . (A6) 



5t 5t v x (5t) 2 v a 



(A7) 



26 



and for a surface of small periodic length a one can approximate 

dx a 



- dt - a- (A8) 

Taking the infinitesimal limit, Eq. (1A6I) becomes 

5E = pEka (^) 2 ^dt. (A9) 

Equation flA9j) identifies the power loss spectrum, which can be evaluated from Eq. flA4j) 

— — = fiEid — j— 2 — — . A10 

dt vl 4cosh 4 (rv t/2) V ' 

Integrating Eq. (lAlOj) over all times along the path of motion and recognizing that the total 



scattering angle for the grazing-angle collision is 9 « 29 i leads to a result that is cubic in 
the incident angle, rather than the quadratic form suggested by Eq. (1A5|^ ^: 

AE = - fiEiTa9f . (All) 
3 

The result of Eq. ( 1A11I) can be understood in physical terms through the following simple 
argument: in a grazing-angle collision the interaction can be assumed to be spread roughly 
equally over sequential and incoherent collisions with A^ lattice sites in a line along the 
surface, and clearly for such grazing-angle conditions A^ oc If each of these individual 

collisions has a fractional energy loss equivalent to Eq. (IA6I) and if the fractional angular 
deflection is the uniform value 59 w (9i + 9f)/N ~ 29i/N then the sum of all energy losses 
over the A^ small collisions along the trajectory is 

AE = . (A12) 

Since A^ is known to be inversely proportional to 0j, Eq. (IA12I) gives the same cubic depen- 
dence on incident angle as the more exact result of Eq. ( 1A11I) . This argument also identifies 
an approximation for estimating the number of surface cells lying along the path of the 
collision:-*^ 

N ~ m: • (A13) 
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2. The lognormal angular distribution 



If the surface potential is exponentially repulsive as in Eq. (fT2j) and additionally if each 
element of the surface is subject to a random displacement in the vertical direction the 
distribution of final scattering angles is lognormal.- 1 ^ This is perhaps most easily understood 
through noting that in a grazing-angle scattering event with V{ Z « v x the final angle is 
approximately given by 



« ~ _J£ 

v~ 



2E z /m 



(A14) 



and this can be related to the classical turning point through 



e 



f 



m 



-Fz /2 = Q- f e -Fzo/2 _ 



(A15) 



Now suppose that the surface is subject to random displacements u z in the vertical direction 
whose probability distribution is Gaussian and given by 



P(u, 



1 



V2 



cxp 



7TO~ z 



u: 



2a 2 



(A16) 



Even though ultimately the displacements u z may be due to thermal motion of the surface 
atoms the motion of the projectile is so fast that they appear static, i.e., the collision time 
is much shorter than a typical phonon period. 

One can then obtain the mean value of Of through 



Of 

which gives the result 



du z e- r ^-«,)/2 



V2 



exp 



TTU 



2a 2 



3 r 2 <r 2 /8 



If the integration variable in Eq. (|A17[) is changed to Of according to 



In 



the result is 



(Of) 



1 

r^7 V vra 2 



dOf exp 



T 2 a 2 



(A17) 



(A18) 



(A19) 



(A20) 
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However, relating Eq. (1A20j) to the standard form in terms of the distribution function P{9f) 
given in the small angle approximation by 



0i 



dO f e f P{9 f ) 



(A21) 



identifies P{9f) as the lognormal distribution 
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(A22) 



For an ordered surface, the roughness will be caused by thermal motion. So even though 
the surface appears static to the rapidly moving projectile, the width of the Gaussian dis- 
tribution can be associated with the mean square thermal displacement 



<«;> • 



(A23) 



APPENDIX B: CONVOLUTION OF MULTIPLE COLLISIONS WITH THE 
SURFACE 



The differential reflection coefficient describing classical scattering is Eq. 



w(k f ,ki) 



h 2 k 



mL 



h 



Ank B T s AE 



exp 



T 7 



(Bl) 



(E f -E t + AE 
4k B T s AE 



under many conditions, including the case of grazing angle scattering at high energies con- 
sidered here, this becomes very nearly Gaussian in form leading to a differential reflection 
coefficient Eq. (IC2I) of the form 

' ,, .,- v ,. " } \ ■ (B2) 



4nk B T s AE 



dE f d 2 VL f 

When this happens the most probable energy is given by the Baule relation for an elastic 
two-body collision Ef = /(//, 9)Ei where \x is the mass ratio, 9 is the scattering angle and 

fM = i v ^ - r; ; 1 i . (B3 ) 
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The recoil energy evaluated at the most probable final energy is AE = <7ta(/-*, 0)Ei with 

g TA {^o) = fi (i + f(fi,e)-2y/fhi,e) cose) . (B4) 

The mean square energy deviation is (Ej) = 2g(/j,,6)EikBTs where 

tfGM) = — • ( B5 ) 



fi cos 9 



The trajectory approximation to g(fi,9) is grAil^, 9) and clearly 

/CM) = 1 - QtaM , (B6) 

and for small scattering angles 

g(v,e) ~ qtaM « , (B7) 

The interaction with the surface is described in terms of successive collisions with N 
surface atoms and each of these collisions has a scattering angle which if averaged over the 
whole path is approximately 9 n = 29i/N since the final scattering intensity is predominantly 
along the circle of energy conservation. The differential reflection coefficient for the n-th 
collision is given by 

d 3 R(k n , k w _x) 1 \ (E n — fnjp, 9 n )E n _i) 2 | 

dE n d?n n ~ ^a n eXP \ al )■ ^ ] 

where the definition of a n is obvious from Eq. (1B2I) . The differential reflection coefficient 
resulting from this series of collisions is the convolution of N different successive collisions 

*m,K> _ f dW j^M... / dEl ^f^M,^M . (B9) 



dE f d 2 n f J dE N d?tt N J 1 1 dE 2 d?tt 2 dExd 2 ^ 

The usual rule that the convolution of Gaussians is a Gaussian whose squared width is 
the sum of the squared widths of each element of the convolution becomes slightly more 
complicated in this case because of the factor f n E n -i in the argument. In this case the final 
result is 

dE f d*n f ~ 0Fcr T exp ] 4 ' ' ,blU) 
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where 



n=l 



h = n /»— > 1 - N n ^) > ( Bn ) 



and 



N N ~ j f29 \ 2 

n=l j=l ^ ' 

where the expressions on the RHS of Eqs. (IBllj) and (1B12I) are obtained after making the 
small angle approximation that all f n are nearly equal to unity and a n are small, and both 
are independent of n. Recalling the definition of AE from Eq. ( 1301) the differential reflection 
coefficient now becomes that of Eq. fl3"TT) : 



APPENDIX C: EXPERIMENTALLY MEASURED QUANTITIES 

The physical quantity measured by these experiments is the differential reflection coeffi- 
cient per unit final energy and solid angle denoted as dR(kf, kj) / dEfdVtf. This is obtained 
from the transition rate by first forming the differential reflection coefficient per unit volume 
of wave vector space 

d 3 R(k f , k<) w(k f ,kj) 

dk f Afc 3 ji ' { } 

where in terms of the quantization length L in each Cartesian direction the element of 
volume in phase space is Afc 3 = (27r/L) 3 and the incident flux normal to the surface is 
ji = hk iz /mL. Changing to spherical coordinates according to dkf = kjrdkfdQf and recalling 



that Ef = h kj/2m leads to 

dR^ki) _ LV w(k/jki) . (C2) 



dE f d 2 VL } {2-nhf k { 
The channel plate detector used in the experiments measures the integral over all final 
particle energies and thus the measured quantity is the differential reflection coefficient per 
unit final polar and azimuthal angles which is expressed as 

d 2 R , n s r d 3 R(k f ,ki) 

where the polar angle used here is measured with respect to the surface plane, and not the 
usual one measured with respect to the surface normal. 
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LIST OF FIGURES 



1 The scattered intensity as a function of 0/ with 6j lying on the Laue circle 
for a 200 eV beam of 3 He incident on LiF(001)< 110 > at an incident angle 
$i = 2° with respect to the surface. All other parameters are specified in the 
text. The broad features are the calculations for the nearly classical case of 
Eq. (150]) and the sharp vertical lines at the diffraction peak positions represent 
the intensities of the quantum diffraction peaks calculated from Eq. (|37|) . 

2 For the same conditions as Fig. [T]the differential reflection coefficient is plotted 
as a function of Of for 0j fixed at the first and third order diffraction peak 
position in the upper and lower panels, respectively. 

3 The energy- dependent differential reflection coefficient of the nearly classical 
model of Eq. (1501) is plotted as a function of final energy Ef at the specular 
and second order diffraction peak positions in the upper and lower panels, 
respectively 
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FIG. 1: The scattered intensity as a function of <f>f with Of lying on the Laue circle for a 200 
eV beam of 3 He incident on LiF(001)< 110 > at an incident angle 0\ = 2° with respect to the 
surface. All other parameters are specified in the text. The broad features are the calculations for 
the nearly classical case of Eq. (|50p and the sharp vertical lines at the diffraction peak positions 
represent the intensities of the quantum diffraction peaks calculated from Eq. (|37p . 
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FIG. 2: For the same conditions as Fig. [T] the differential reflection coefficient is plotted as a 
function of Of for (j)f fixed at the first and third order diffraction peak position in the upper and 
lower panels, respectively. 



36 




199.95 



200.05 



E f (eV) 



FIG. 3: The energy-dependent differential reflection coefficient of the nearly classical model of 
Eq. (|50p is plotted as a function of final energy Ef at the specular and second order diffraction 
peak positions in the upper and lower panels, respectively. 
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